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Abstract 

We study the temperature dependence of the rj and rj 1 meson masses within the frame- 
work of U(3)l x U(3)r chiral perturbation theory, up to next-to-leading order in a simul- 
taneous expansion in momenta, quark masses and number of colours. We find that both 
masses decrease at low temperatures, but only very slightly. We analyze higher order 
corrections and argue that large N c suggests a discontinuous drop of M^t at the critical 
temperature of deconfinement T c , consistent with a first order transition to a phase with 
approximate U{1)a symmetry. 
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1 Introduction 



The fate of the U{1)a symmetry of QCD at finite temperature is a fascinating problem |TJ, fj. 
which could also have interesting consequences for the ongoing heavy ion collisions program 
and, possibly, for cosmology @, ||, || [7], ||. Even a partial restoration of U(1)a symmetry in 
the vicinity of the critical temperature of chiral symmetry breaking (T c ~ 200 MeV) could 
dramatically change the mass and mixing pattern of the lightest neutral mesons (tt°, t] and 77'), 
with signals including enhanced strangeness production or the more speculative possibility of 
forming parity violating disoriented rf condensates in heavy ion collisions |7|, |9|, |10|] . 

Our aim in the present paper is rather modest: we will study the shift of the mass of the t] 
and rj' mesons at low temperatures, in a regime in which the hadronic gas is mostly composed 
of pions. We will work in the framework of U(3)l x U (3)# chiral perturbation theory (xPT), in 
a simultaneous expansion in momenta, quark masses, number of colours N c , and temperature 
T. 

Our motivation for doing this investigation was threefold. First, the predictions of xPT in a 
pion thermal bath, although limited in scope to T < few are essentially model independent 
(see for instance, the review of Smilga Jll|] and references therein). Given the phenomenological 
success of the large N c expansion in vacuum, one might perhaps hope that the predictions of 
the present work are as robust. Next, we wanted to see to which extent the results derived in 



Ref. |T2[ could be amended. As us, the authors have computed the shift of M v and M r] i at low 



temperatures using the Di Vecchia-Veneziano-Witten effective lagrangian (DVW) ]T3|, II], |T5 
but only to leading order AM 2 , ~ T 2 . However, it was not clear to us whether the leading order 
DVW lagrangian was a good approximation for this problem. Although the parameters of the 
lagrangian can be fitted to the observed mass and mixing pattern of the r\ and rj' mesons to 
within 10% [16|, |Tj]], the decay rates predicted for rf — > r)7nr are off the experimental values by a 
factor of about 40. This issue, which is obviously relevant in order to determine the shift of M v > 



in a pion bath, is however easily cured at next-to-leading order in the large N c expansion [18|, [19] 



As we will show, next-to-leading order corrections are also quite important at finite temperature, 
but not to the point of dramatically changing the conclusion of Ref. [12]: at low temperatures, 
Mrf stays essentially constant. Finally, we wanted to see what the large N c expansion could 
teach us about the fate of the U(1)a symmetry at finite temperature. At zero temperature, in 
the confined phase of QCD, large N c arguments predict that Mi oc l/N c . On the other hand, 
at very high temperatures, in the quark-gluon plasma phase of QCD (T ^> /ihadr ~ 200 MeV), 
instanton calculus is reliable and predicts an effective restoration of U(1)a symmetry. Because 
of screening, instanton effects are suppressed at very large temperatures oc exp(— 8n 2 /g(T) 2 ). 
At large N c , the suppression is more important, as 1/g 2 — > N c /X with fixed 't Hooft coupling 
A = g 2 N c , and the exponential tends to vanish exp(— N c /X) — > as N c increases. Because of 
asymptotic freedom, A growths at lower temperatures and the instanton argument breaksdown. 
However, for N c large enough, a natural assumption is that the exponential suppression holds 
all the way down to the critical temperature of deconfinement T c ~ /^hadr [0- Although we 
have no proof of this statement, such a behaviour seems natural given the large release of 
entropy cx N 2 at T c and is actually known to occur in models in two dimensions |2T) |. With 
this assumption, M v i can be taken as an order parameter for U(1)a symmetry restoration at 
T c . In Ref. 0, some information on the behaviour of M n i near T c could be extracted assuming 
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that the deconfining phase transition could be of second order at large N c We will argue 
here that large N c favours a sharp drop of M v > at T c , consistent with first order transition to 
the phase with (approximate) U{1)a symmetry. 

Our paper is organized as follows. In the next section, we briefly review U{3)l x U(3)r 
xPT, which extends the framework of the large N c DVW effective lagrangian beyond leading 
order. For definitiveness, we refer to the recent analysis of Herrera-Siklody et al. [22j. We then 
discuss the implications of these corrections at low temperature, in presence of a pion thermal 
bath. Most formulas are relegated to the appendix. In the last section, we speculate on the 
effect of higher order corrections in the large N c expansion and draw the conclusions. 



2 Sketch of U(3)l x U(3)r chiral perturbation theory 



At low energies and temperatures, the dynamics of QCD is governed by an approximate 
SU(3)l x SU(3)r chiral symmetry which is spontaneously broken to the diagonal SU(3) in 
vacuum. If the mass of the up, down and strange quarks were vanishing, the symmetry 
would become exact and there would be eight massless Goldstone bosons. Phenomenologi- 
cal lagrangians, which treat the mass of the quarks as small perturbations, provide a powerful 
framework, known as Chiral Perturbation Theory (%PT), to study the properties of the light 
7r, K, and 77 mesons p3| , The rf meson doesn't a priori fit in this frame. It is substantially 
heavier than the other eight light mesons, and, in vacuum, would stay so even in the chiral limit 
of zero quark masses, because it receives most of its mass from the U(1)a anomaly through 
non-perturbative instanton-like effects [25|. The effect of the axial anomaly can however be con- 
veniently turned off by going to the limit of large number of colours N c |27|, ^] . At infinite 
N c and in the chiral limit, the global symmetry becomes U(3)l x U(3)r, spontaneously broken 
in vacuum to U(3), with nine massless Goldstone bosons. Like chiral symmetry breaking effects 
by finite quark masses, 1/N C suppressed contributions can be systematically introduced as per- 



turbations in an effective lagrangian, an approach which has been quite fruitful [fT3|, \L% |19[ |29 [. 
A systematic analysis of next-to-leading corrections, including 0(p A ) operators, has been initi- 
ated in the recent 0, 0. We refer to these latter works for more details and follow their 
conventions for ease of reference. We will work in Euclidean spacetime with metric g^ u = 
and use the imaginary time formalism to compute the thermal corrections. 



2.1 Leading order 



The leading order effective lagrangian is well-known JT3], [14], In the notation of it is 
written as 



C 



LO 



where U is the U(3) matrix 



(v 02 X 2 + (dprfdpU) - (U^ X + x ] u) 
U = e^/ f , 



(1) 
(2) 
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with $ the pseudoscalar meson matrix 



(3) 



and / = f n = 92.4 MeV at leading order. The mass matrix is 

X = 2 B diag(m u , m d , m s ) , (4) 

but we shall neglect isospin breaking effects (m u = rrid = m). The constant B is related to the 
value of the (qq) condensate, M 2 = 2mB = —2m(qq)/f 2 at leading order. The combination 

X(x) = (log U(x)) + i9qcd = i —Vo + ^qcd , (5) 

is invariant under f/(3)iX U(3)r transformations, (log£7) — > (log^ij?/^)) + 2i(a) and #qcd — ► 
$qcd — 2(a). Because of this, any arbitrary function of X can a priorienter in the construction of 
the effective lagrangian, with thus little predictive power. This is where the large N c expansion 
comes to the rescue by limiting the number of operators that can contribute at each level of 
approximation. In the chiral limit (m,m s — > 0), Eq. ([I]) gives 

Mj = -3v 02 , (6) 

which is the celebrated Veneziano-Witten relation for three massless flavours p7| , p8fl , with 
vq2 = —^r/f 2 ~ l/iV c , where r is the topological susceptibility of pure Yang-Mills theory. The 
rationale of U(3)l x U(3)r xPT is to count powers of p 2 , m q , and l/N c on the same level 

0(6) ~p 2 ~ ~ 1/N C . (7) 

According to this counting rule, the leading order lagrangian (|l|) is O(5 ) because f 2 ~ 0(N C ) []. 

At leading order, there are four unknown parameters in the lagrangian: /, t>02, an d the 
combinations mB and m s B (or x = m s /m — 1). On the other hand, we have at our disposal 
seven observables: /„-, Jk, the four masses of the light mesons, and the r]-r]' mixing angle 9. 
Using M v r as input and the formula? given for reference in the appendix, one obtains 

f2 _ f2 _ f2 

J — Jtx — J K > 

2mB = M 2 , 



which predict that 9 ~ —20° and 



x ~ 24.1 , 

v 02 ~ -0.22 GeV 2 , 

M v ~ 494.4 MeV . (9) 



1 Note that the field expansion of U brings further powers of 1// ~ l/yN c . The 0(5) counting is to be 
understood to hold at the operator level. 
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Remarkably, the latter number is only 10% off the experimental value M v = 547.3 MeV. It 
is however known that adjusting the parameters cannot improve the prediction because the 
ratio M^/M^, has an upper boundf] |l(|. One has to take into account next-to- leading order 
corrections to reach agreement fl~7|j . 

At leading order, the only coupling between rj and the pions is from the quark mass term 
in the lagrangian (|l|) and is thus chirally suppressed. The amplitude for 7/ — > t]tiit is then 



A = ^(2v / 2cos(20) - sin(20)) . 



The corresponding decay rates 

r(r/ -> T]n°n°) = 1.0 keV , 

T(rf r]7v + 7i-) = 1.9 keV « 2 x T(r]' 

are however much smaller than the experimental ones, 



r]7i°7i°) 



r exp (V 



42.0 ± 4.2 keV 
88.9 ± 7.6 keV 



(12) 



(13) 



(14) 



We will not speculate on the reasons for this well-known discrepancy (see Refs. [T9|, |31] for a 
more recent discussion), but simply note that within the present framework, this issue can also 
be resolved at next-to-leading order [TR [T9T]p|. 



2 Assuming M v — to simplify, Eq. (Q) gives 



Ml = 3 - ?7 - y/9 + 2y + y 2 
M 2 3 - y + V9 + 2y + ' 

where y = 9vo2/2(Mj c — M%). This ratio reaches a maximum at y = — 3 (note that W02 < 0) corresponding to 

^< 0.518, (H, 

to be compared with the measured ratio M^/M^i ~ 0.571. Taking into account M T ^ improves things, but 
not enough. 

3 Note that the amplitude ( p^ ) is constant and vanishes in the limit m u — md = 0, for any m s . However, 
general arguments p9| (and a fit to experimental data) indicate that for m u — md — 0, the amplitude should 
behave like A — const x ■ , where pi are the momenta of the outgoing pions, and where the constant 
is vanishing as the strange quark mass goes to zero. As shown in Ref. |l9| ], this behaviour can be easily 
accommodated by introducing higher-order terms, an approach that is systematized by the 5 expansion fl8| . 
The smallness of the leading order contribution is then considered as a mere accident, related to the smallness 
of the ratio M%/M^. 
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2.2 Next-to-leading order 

In our case, at next-to-leading order, 0(5), only a few more terms can be added to the lagrangian 

+ 4 ( - v 31 X(W X ~ X ] U) + v i0 (Wd,U)(Wd,U) 

x (15) 

+ iv 50 (U 9fj,U) d^QCB + V 60 ^^qcd^^qcd) 



- M O -M 3 O 3 + L 5 O 5 -L 8 O 8 , 
where the 00,3,5,8 are 0(p 4 ) operators whose coupling constants are 0(N C ): 

O = (dpUduWdpUdvtf) , 

5 = (d,Wd,u(W x + x j u)) , (16) 
o 8 = (x ] u x ] u + W x rfx) ■ 

The couplings f 40 , v 50 , and v 60 are not independent and either one of them can be set to 
zero by an appropriate change of variables, rj / f — > rj /f + k9qcb- We shall choose t> 40 = 0. 
Moreover, f 50 and Veo will not appear in our calculations and can be discarded. At 0(S), the 
only coupling related to the breaking of U(1)a symmetry is thus V31 ~ 0(1/N C ). Note that the 
corresponding operator is also chirally suppressed oc m q . 

At next-to-leading order, seven unknown parameters enter in the definition of the meson 
mass matrix: /, t> 2, ^31) ^5,8, together with the quark masses m and m s (see appendix for 
details). These can be expressed in terms of seven independent observables: /„-, fx, M n , Mk, 
M v , ML and the rj-rj' mixing angle 9 At this level, large A^ c xPT is thus not predictive. 
The strategy adopted in Ref. |30| was to impose that 0(8) corrections are not too large so 



that the large N c expansion makes sense. For mixing angle in the range 20° < 6 < 24°, the fit 
gives |30[ 

0.980 < ImBjMl < 0.988 , 
18.3 < x < 20.9 , 
0.214 GeV 2 < |voa| < 0.239 GeV 2 , (17) 
1.35 10~ 3 < U < 1.57 10- 3 



-0.164 < u 31 < -0.161 . 

together with / = 90.8 MeV and L 5 = 2.0 10~ 3 which are fixed by f n and fx- Note that if v 02 
does only change by about 10%, the shift in m s is quite large, ~ 20-25. 

Because they have four derivatives, the operators Oq and O3 do not contribute to the 
meson mass matrix in vacuum. However, they give the dominant contributions to the decay 
7/ — > i]7i7i ||18|| . This is essentially because the extra derivatives introduce large amplification 
factors, oc (M v > / M n ) 2 , with respect to the leading order amplitude^. The observed decay rates 

4 This may actually casts some doubts on the validity of xPT for such processes as one could expect higher- 
order effects to give non-negligible contributions to the decay r/' — > rjirir. One may nevertheless hope that the 
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are well reproduced with 

M ~ 1.2 10" 3 , 

(18) 

M 3 ~ -0.4 10^ 3 , 

values which can be independently inferred from the known L±, L 2 , and L 3 of SU(3) x SU(3) 
xPT (in the nomenclature of Gasser and Leutwyler [P^])p|. 

Thus all the parameters of the next-to-leading order effective lagrangian are fixed by low- 
energy phenomenology. 



3 M' in a pion thermal bath 



In Ref. [0, the leading order lagrangian ([]]) has been used to study the shift of M v and M v / 
at one-loop in a pion thermal bath. The effect they found is very tiny, as at T ~ 200 MeV the 
relative mass shifts are only 0(0.1%). The reason for this is easy to understand. The r\ and r{ 
mesons receive most of their mass from the topological susceptibility term oc t>02 and/or from 
the strange quark mass, while the pion thermal corrections only modifies the tiny contribution 
from the pion mass term cx M%. Thermal kaons could give a larger effect, cx Mj^, but the 
density of these is exponentially suppressed at low temperatures, oc exp(— Mk/T). One might 
wonder whether next-to-leading order corrections could directly affect the contribution of the 
leading order U(1)a breaking term vq2- As we have seen in the previous section, five extra 
operators appear at next-to-leading order in the large N c expansion and, of these, only the 
one with coupling t> 31 is related to U(1)a symmetry breaking. Unfortunately, this term is also 
chirally suppressed, cx m, and its contribution is only 0(v 3 iM 2 T 2 / f 2 ). At temperatures of 
interest, this is small compared to Vq2, but of the same magnitude as the leading order thermal 
correction. The other four operators will also contribute, but in a less interesting way, as they 
are invariant under U(1)a- Furthermore, their effects are also oc M 2 . 

We have computed the shift of the mass of rj and rf at one-loop, at next-to-leading order in 
the expansion in 5. We have not taken into account two-loop corrections from the leading order 
lagrangian. Although it is not clear whether this is legitimate numerically speaking, neglecting 
these is however consistent with the rules of large N c chiral perturbation theory. Indeed, the 
natural extension of 5 power-counting to finite temperature is 

0{5) ~ p 2 ~ m q ~ 1/N C ~ T 2 . (19) 

At leading order, M 2 , = 0(6) ~ 1/N C and the one-loop thermal correction is oc M 2 T 2 /f 2 ~ <5 3 . 
At two-loop, using the leading order lagrangian, the shift is oc M 2 T 4 / f 4 ~ <5 5 , while at one-loop 
using the next-to-leading order lagrangian, the shift is typically v 3 iM 2 T 2 / f 2 ~ 5 A (using v 3 i ~ 
l/iV c ) and thus dominant. Consistency thus requires to neglect the two-loop contributions. This 
greatly simplifies the calculations which are a bit cumbersome, but otherwise straightforward. 

large N c expansion is still reliable and that these corrections are 1/N C suppressed. Whether this is true is 
unfortunately hard to check as we would evidently have too few hadronic data to completely fit the parameters 
of the effective lagrangian at higher orders in the S expansion. Of course, this is precisely why the large N c 
expansion is invoked in the first place. 

5 According to Ref. @, M = §(Li + L 2 ) + O(iV c ) and M 3 = L 3 + 2M . 
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Figure 1: One- loop pion corrections to M v and M v / at low temperature. The black box represent 
insertions of next-to-leading order operators. 
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Figure 2: Leading order (solid line) and next-to-leading order (dashed line) contributions to 
M V (T). Both curves are normalized to M v (0). 



The relevant diagrams are those of Fig. |l|, where the loops contain only pions. At next-to- 
leading order there are two related thermal loops: 

/i(r) ^jl/(S™' (2o) 

with K 2 = + k 2 , where ko = 2irnT, with n integer, are the Matsubara frequencies and 

00 r rl 3 k K 2 

(21) 

with hiT) = —M 2 Ii(T). As usual, we drop the ultraviolet divergent part of the pion loops as 
these can in principle be reabsorbed in vacuum parameters, including next-to-next-to-leading 
order counter-terms. The sum over n can then be readily evaluated using standard tech- 
niques p2f , 

h[T) =t^—^^±l K A, (22) 

J (2ny uj exp(u)/T) — 1 2n 2 ^ n v T ' 
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Figure 3: Leading order (solid line) and next-to-leading order (dashed line) contributions to 
M V >(T). Both curves are normalized to M v >(0). 

where uj 2 = k 2 + M 2 . For instance, for T > M n 

h(T) « g . (23) 

In the sequel we simply compute fl22|) numerically. As the relevant formula are not particularly 
transparent, we have relegated them to the appendix. Fig. H and Fig. [3] show the shift of 
and M v at low temperature both at leading and next-to-leading order. To be definite we 
have chosen the set of parameters corresponding to 9 ~ —20°. The net thermal effects are 
not dramatic: both masses decrease, but only slightly. As expected, the shift of the mass 
M ! n is more pronounced at next-to- leading order, but the effect is not very significant. Again, 
this is because, both at leading and next-to-leading orders, the thermal corrections are chirally 
suppressed, oc M 2 T 2 / f 2 . For completeness, we have also plotted in Fig. f| the shift of the 
mixing angle at low temperature. As both the rj and rf masses diminish, the angle is not very 
much affected. It decreases a bit (toward ideal mixing?), consistent with the relatively larger 
shift of M v >. 

4 Lessons from large N c ? 

As we have seen in the previous section, the mass of rj' is almost not affected at low temper- 
atures in a pion bath. This is because, at this order, the pion thermal corrections are chirally 
suppressed, smaller than M 2 ~ 0.02 GeV 2 , and thus essentially negligible compared to the 
contribution from the U(1)a symmetry breaking term v 02 — 0.22 GeV 2 . In particular, in the 
chiral limit, M n = 0, all the corrections vanish and M v > is temperature independent up to 
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Figure 4: Leading order (solid line) and next-to-leading order (dashed line) contributions to 
tan2#(T). Both curves are normalized to tan 26(0). 



next-to-leading order in xPT. In the chiral limit, the leading contribution from pions to the 
shift of M v > presumably arise from 0(p A ) operators like[] 

C ~ -J- X 2 (duUdvrfdpUdvrf) . (25) 

The coupling is 0(1/N C ) because there is a factor of 1/N 2 coming with X 2 and one of N c 
from the coupling of the 0{j) A ) operator. The best way to see this is to replace the coupling 
M ~ iV c (or M 3 ) from the next-to-leading order lagrangian by a function of X, M — > M (X) 
and expand to second order in X, which brings down a factor of 1/N 2 . Because there are at 
least four pions in the expansion of the operator (P5|), the leading pion thermal correction to 
M v i in the chiral limit is a two-loop effect, 

1 T 8 IT 8 
SM 2 (T) ~ — — r^—— — (26) 

We have made the large N c dependence of the pion decay constant f 2 manifest by defining 
f 2 ~ A^c/ihadr- Of course, the sign of the correction is not known and M v i could go up or 

6 The 0(p 2 ) operator 

£~f 2 ^(d M U%U) , (24) 

can contribute at one- loop if and only if M% ^ 0. It could contribute at two-loop order in the chiral limit oc T 4 , 
but doesn't because pion interactions are too soft. This is a well-known feature of pion thermal corrections 
which is for instance manifest in the absence of T 6 terms in the free energy of a pion gas in the chiral limit J33| , 
or in the fact that massless thermal pions move at the speed of light to order T 2 Q] . 
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down. Also, if we compare with MB(0) ~ /^hadr/^c, we infer that M v i is quasi-constant for 
temperatures 

T « % ~ iVfVhadr • (27) 

In the large iV c framework, the natural scale for deconfinement is T c ~ /ihadr, which is also 
the temperature at which the pions from the hadronic gas overlap. It is natural to assume 
that chiral symmetry restoration takes place at the same temperature, driven by the release of 
0(N 2 ) gluon degrees of freedom p], |35|. The estimate in Eq. ( 2% then seems to imply that 
is essentially constant up to the temperature of deconfinement, since T* 3> T c ~ /ihadr for N c 
large. This conclusion is however premature because the low momentum expansion breaks down 
near T c and we must take into account the contribution of operators with arbitrary number of 
derivatives. We claim that the dominant operators at large N c are of the form 

£ ~ TtAct^ 2 {dnUdnrf . . . d^Ud^d^U . . . dpjrf) . (28) 

iV c /^hadr 

These operators are irrelevant at low energies but become marginal for d ~ T c . A six-derivative 
operator, for instance, first contribute at three loops <5M 2 , oc T 12 /(A^/i h ° dr ). For comparison, 
the contribution of a three-loop diagram with a four-derivatives (NLO) and a two- derivatives 
vertices (LO) is oc 1/N C T 10 / f 8 ~ l/A^T 10 //i^ adr and is subdominant at large N c . For generic 
k, the operators of Eq. (p8|) give £M 2 , ~ T 4fc /(^ +2 / u^ dr 2 ). The ratio of two consecutive terms 
k and k + 1 becomes 0(1) at T ~ iV^Vhadp, independent of k. At large N c , T c < T' < T* 
and the perturbative expansion still breaks down above the temperature of deconfinement. 

Another set of operators could be relevant at large N c because the rf is then rather light, 
M v ' ~ /^hadr/-^ r c 1/ ' 2 < T c . Thus one should include operators that involve arbitrary powers of the 
f] 1 field, like 

C ~ PF (§-) (d,U%U) 

/ 2 4 \ l^yj 

= (wk + ^fe + • • •) (<W) 2 + P ion terms > 

which contributes to the wave-function renormalization of 7/, 



5Z V , ~ T 2 /(N^U) + T'/iN^U) + . . . , (30) 



or terms of the form 



£ ~ N^UG ( ~ 4^% 2 + #i4 % 4 + #27T^ % 6 + • • • • (31) 



NJ N c ~ ~ "^Ldr^c 7 

However, a common feature of these operators is that they are very suppressed at large N c . 
They become important only for T ~ N^ 2 ^^, much higher than T' so that their contribution 
is subleading compared to operators like in Eq. (ffif). 

Can we conclude anything from these considerations? In all the cases discussed above, the 
leading thermal corrections to M v ', in the chiral limit and for N c large, become important 
for temperature which are higher than the critical temperature of deconfinement T c ~ /ihadr 
by a factor of iVJ. Although the value of 7 is hard to guess, as various corrections can get 
mixed up, we believe it is reasonable to conjecture that 7 is strictly positive. This implies that 
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just below T c , M V >(T) = M^(0) to a very good approximation. The standard lore is that the 
deconfining phase transition at T c is of first order for N c large |$3[} Because the temperature 
at which hadronic interactions can affect M v > is (very much) larger than the temperature of 
deconfinement , we expect that changes in M v > will be instead triggered by the release of the 
large number of gluons and will thus drop dis continuously at T c , i.e. that there is a first order 
transition to a phase with (approximate) U(1)a symmetry. 

This behaviour is not inconsistent with various other expectations. For three light quark 
flavours, Nf > 3, the transition to the chirally symmetric phase is probably first order while 
for Nf = 2, the phase transition is supposed to be of second order, in the universality class of 
0(4) J2[. It has been argued by Smilga that the latter behaviour is not inconsistent with a first 



order deconfining phase transition at large N c [IT|. The reason is that, unlike for M v >, there is 
an infinite subset of thermal corrections that contribute to the same order in N c to the shift of 
the quark condensate S = (qq), 



S(T) =£11- rtj^^F . (32) 

Even though thermal corrections are suppressed like the (unknown) function F(x) may 

be singular near, but below T c ~ yUhadr- If, for instance, F has a simple pole at T = T < T c , 
F ~ /^hadr/ — ) an< i the chiral phase transition is second order with a critical region near 
T c , that is of order 

Alternatively, if U(1)a symmetry is effectively restored at T c , the large N c behaviour (^) is 
also consistent with a fluctuation induced first order phase transition. For Nf = 1 finally, chiral 
symmetry is broken by the anomaly at all temperature and there is no chiral phase transition. 
However, if instanton transitions are strongly suppressed just above T c , chiral symmetry can 
be effectively restored and the phase transition is presumably first order. 



5 Conclusions 

We have studied the behaviour of the mass of the rf pseudoscalar meson at finite temperature 
using constraints from chiral symmetry and large N c power counting. The main conclusion to 
be drawn from this work is that is essentially unchanged at low temperatures. A tentative 
analysis of the effect of leading higher order corrections at large iV c suggests that M n * changes 
discontinuously at the temperature of deconfinement. The implications of these considerations 
for the real world, i.e. N c = 3, are not quite clear as we would expect the suppression of 
instanton effects only at asymptotically high temperatures. It is however striking that recent 
lattice simulations, with Nf = 2 staggered [B^] and domain wall [35] fermions, both show a 



strong suppression of U(1)a breaking effects at low temperatures T ~ 1.2 T c . Because this 

7 Various arguments, including recent developments in string theory (see Sect. 6.2.2 in Ref (36|) and lattice 
simulations of N c = 4 pure Yang-Mills theory |j7f , favour a first order deconfining phase transition. A case for 
a second order phase transition has been made in Ref. [21j], in light of the structure of the Columbia diagram. 
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Nc>3 



Tc 



Figure 5: Plot of the 7r and 5 (ao) near the critical temperature for N c = 3 (continuous lines) 
and possible changes for N c ^> 3 (dashed lines). 



temperature is outside the critical region, the order of the chiral phase transition is probably 
not affected. It could be of interest to consider doing simulations with N c > 3, although this 
would probably be time consuming, or maybe with one flavour and various N c . Consider for 
example a plot of the tc and 5 susceptibilities near the critical temperature as computed on 
the lattice IPBfl . Large N c arguments suggest that the curves of the susceptibilities would be 



flatter below T c — because the confined phase is colder — and that splitting between n and 5 
(a.k.a. a ) should be narrower above T c — because U(1)a breaking is more suppressed, maybe 
like in Fig. [5] — . The Columbia diagram could change accordingly: the critical line around the 
region of small M u = Ma masses would move as in Fig. |6|. 
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Appendix 

5.1 Useful formulae 
5.1.1 Leading order 

2mB = Ml , (34) 
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Nc = 3 



Nc»3 



Nf=0 Nf=2 





Nf=l Nf=3 



Figure 6: Columbia phase diagram as function of M u = Md (horizontal axis) and M s (vertical 
axes) for N c = 3 (left, the blue dot is where QCD stands) and how it could evolve at large N c 
(right). 



x 



Ml 
2—^-2 



Ml 

/ J ft J 

= Ml-^^^Ml-Ml)^ 



(35) 
(36) 
(37) 



5.1.2 Next-to-leading order 

Some definitions: Am, Ajv are defined as 



A M = -(M 2 K -M 2 W )(2L 8 -L 5j 



L 5 

3v 31 - 12— f 02 • 



(38) 
(39) 



Next-to-leading order parameters: The next-to-leading order parameters can be ex- 
pressed in terms of observables through 



2mB 



Ml (l - "L A 



M 2 K -Ml 



M 



M 



x 



K , 



Ml 



1 - A 



f = fll-A^M 2 j , 



(40) 
(41) 
(42) 
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3v 02 = M, 



■2. Ml- + Ml 



+ 



2\/2 

-^-(M| - M 2 )(l + A M - Ajv) tanfl 



{M 2 K - Ml) A M - (2M| + M^Ajv 
1 



/ 2 



4(Mi-M2)' 



where 



M 2 + 3M 2 - AMI + 3(M 2 - M 2 ) sin 2 9 
Am ~ 77172 



A P = ^-1. 



4(M 2 - Ml) 

3 (M 2 -M 2 )sin2fl 
4^2 M 2 K -Ml ' 



A N = 1 + * + Am 



5.2 Results at T = 
5.2.1 Leading order 
Mass matrix: 



m 2 8 = 1 -{AM 2 K -Ml), 

™so = -^(M 2 K -Ml), 

m 2 = ^(2M 2 K + Ml)-3v 02 

= \(2Ml + Ml)- 2 -y{Ml-Ml), 



where 



Mixing angle: 



9v, 



02 



W - *f, 2 ) 



tmM5 = M 

m 00 ~ m 88 1 +V 
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Physical masses: 



Ml - Ml 



Ml = M 2 K - K o "'* ( y -J9 + 2y + y*) , 



with 



Ml + M% = m 2 88 + m 2 m = 2M 2 K - \y(M 2 K - Ml) . 



5.2.2 Next-to-leading order 
Mass matrix: 



^88 = ^Ml-Ml) + -(M 2 K -M 2 )A M , 

m 2 80 = -^f(M 2 K -M 2 )(l + A M -A N ), 

m 2 00 = ^(2M 2 K + Ml)-3v 02 

+ ^(M 2 K -M 2 )A M -^(2M 2 K + M 2 )A N 

= ^(2M 2 K + M 2 )(1-2A N ) + ^(M 2 K -M 2 )A M 

- \y{M 2 K -Ml). 



Mixing angle: 



nn 2^2 / y / 1 2M 2 K + Ml\ N 



Physical masses: 



m 2 v = M 2 K - Ml 3 M h y + ^ + 2y + y 2 ) 

+ (1 9 + y ) (Ml -M 2 ) A M 

1 ( 0M 2 , M 2 S(2M 2 K - 3M 2 ) - y(2M 2 K + M 2 ) \ 
- -[2M K + M„ -j==== )A N 

Ml -Ml 



Ml = M 2 K - K - -* ( y -J9 + 2y + y*) 
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with 



M 2 + M 2 = 2M 2 K -\y{M 2 K -Ml) 



+ 2(M 2 K -M 2 )A M -^(2M 2 K + M 2 )A N . 



5.3 Results at T ^ 
5.3.1 Leading order 
Mass matrix: 



mls(T) = ml 8 (0) - ^-I(T) , 



M 2 

mUT) = ml (0) j=?—I(T) , 



M 2 

mUT) = m 2 00 (0) - -f/(T) , 



V 2 
Ml 
V2fZ 
Ml 
ft 



where 



Mixing angle: 



2tt 2 ^ n 1 V T / ~ ' 12 ' 



t a„ 2 «(T)^an2«(0) + ^^?-^l/(T) 



Physical masses: 

M 



*(T) = M 2 (0) - 3 -^I(T) ( 1 ■ 9 + " 



4/2 v 7 V 3V9 + 2y + yV ' 



MJ(T) = JWJ(O) - -Ml{T) { 1 - ° + U 



4/2 v 7 V 3V9 + 2y + yV ' 
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with 



3 M 2 

M 2 (T) + M 2 (T) = M 2 (0) + Ml (0) - --^/(T) 



5.3.2 Next-to-leading order 
Mass matrix: 



/if 2 / oM 2 f 3 

« r > = <^-wi IiT) [ 1+ w^m( Ap+ 2 AM 

Ml 



2/2 v ^ M£-M 2 

YM 

f 2 



+ 24^f(M + M 3 



- A„ + ,f,M. + M I ,( 1 -^)), 

M 2 / 9M 2 / 3 \ 

m^T) = m§o(0) - -jf-KT) [l + _' M{ (a p + -A M ) 

- 2Aw + 24 |(M + M 3 )( 1 -|M)). 

Mixing angle: 



2M 2 / 3 A \ 2 , 

l-y A 3 2M|-(l + j/ )M 2 

1 + 2/ ^ (l + y)y M|-M 2 1 



M 2 . J Ml - \(M 2 K - Ml 

+ 24-f (M + M 3 ) ' " 5 " 



f2 v u a M M 2 



Physical masses: 



M,?(T) = MJ(0) - |f /(T) (l + j^j^, 



Ml -Ml \\ 3\f¥T2yTW , 
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+ 3 / + 27(3 + y) + (9 + y)y 



2\ 3(9 + 2y + y 2 ) 3 / 2 

\y 2 AM 2 K -Ml" 
"3(9 + 2t/ + ?/ 2 ) 3 / 2 M 2 



A 



4 A (3 + y)((9 + y) + (3 + y)y) 
3 ^ (9 + 2y + y 2 ) 3 / 2 

l % M 2 \ 

^(9 + 2 2 / + |/ 2 )3/ 2 M 2 -M 2 yl " 

+ 24^(M + M 3 ) f l + _-£±^ == 
/. 2V V 3 V9 + 2y + 2/ 2 



9" V v /9 + 2y + y' 2 J M 2 
Mj(T) = Mj(0) - ^/(T) I 1 - 



V ' K ' 4 / 2 1 J l 3 V9 + 2y + y 2 



Ml-M^VV 3V9 + 2y + y 2 / 
+ 3 / 27(3 + y) + (9 + y)y 2 



2\ 3(9 + 2y + y 2 ) 3 / 2 

\y 2 m\ -Ml " 

'3(9 + 2y + y 2 ) 3 / 2 M 2 



A 



4 / (3 + y)((9 + y) + (3 + y)y) 
3 ^ (9 + 2y + y 2 ) 3 / 2 

6y M 2 \ 

(9 + 2?/ + 2 / 2 ) 3 / 2 M 2 : -M 2 y ' N 

•\/:, , 9 + y 



+ 24-f(M + M 3 ) 1 



/r 'V 3 V9 + 2y + y* 

A v (i _ 3 + ^ ^ 
9^ V9 + 2y + y 2 y ' M 2 j 



3 M 2 

JWJ(T)+JWJ(T) = M 2 (0) + M 2 (0)--jfI(T)x 



2Ml ( . 3 . \ 4 



(1 + Afrfl! A -+2 A "J-3 A " 
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+ 24 | ( M + M 3 ,( 1 -i,M_M)). (77) 
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